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Abstract 

The Roberge- Weiss (RW) phase transition in the imaginary chemical potential region is analyzed 
by the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model. In the RW phase transition, 
the charge-conjugation symmetry is spontaneously broken, while the extended Z3 symmetry (the 
RW periodicity) is preserved. The RW transition is of second order at the endpoint. At the zero 
chemical potential, a crossover deconfinement transition appears as a remnant of the second-order 
RW phase transition at the endpoint, while the charge-conjugation symmetry is always preserved. 
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I. INTRODUCTION 



One of the most fascinating and essential subjects in hadron physics is to explore the 
phase diagram of quantum chromodynamics (QCD). QCD is a remarkable theory in the 
sense that it is renormalizable and parameter free. The thermodynamics of QCD is well 
defined, nevertheless not clearly understood because of the nonperturbative nature. A pow- 
erful method of exploring the phase diagram is lattice QCD (LQCD) as the first-principle 



calculation, but it has the wel 



known sign problem when the quark chemical potential (/i) 



is real; for example, see Ref. and references therein. Although several approaches such as 
the reweighting method the Taylor expansion method js] and the analytic continuation 
to the real chemical potential (/xr) from the imaginary chemical potential {nij 4j-|ll| have 
been proposed, these are still far from perfection. 

So far the phase diagram in the /xr re gion has been analyzed by effective models such 
as the Nambu-Jona-Lasinio (NJL) model 12h20|| and the Polyakov-loop extended Nambu- 
Jona-Lasinio (PNJL) model 



2J4I 



41|. The NJL model describes the chiral symmetry break- 



ing, but not the confinement mechanism. The PNJL model is designed [2j] to make it 
possible to treat both the mechanisms. It is reported that the confinement mechanism shifts 
the critical endpoint of the chiral phase transition toward larger T and smaller /iR 29.I35I.I37I. 
When the chemical potential is imaginary, that is /i = ifij = iT6, LQCD has no sign 



problem, so that LQCD data are available there y|-lll|. The first essential work on the fxj 
region was made by Roberge and Weiss (RW) 42]. They found that the thermodynamic 
potential i^qcol^') has a periodicity, i^qcol^') = ^qcd{(^ + 27rA;/3), for any integer k. The 
RW periodicity was proven by showing that f^qcoi^ + ^nk/S) is reduced to ^?qcd(^) with 
the Z3 transformation. 



q Uq, UA^U-^ - i/g{dM)U 



r-l 



(1) 



where q is the quark field, A^, is the gauge field and U{x, r) are elements of SU(3) with 
U{x, 1/T) = exp{—2ink/3)U{x, 0). This means that i^QCD(^) is invariant under the extended 
Z3 transformation jsG], 



9^9 + 2TTk/3, 

q Uq, A^ UA^U-^ - i/g{dM)U-^ . 



(2) 



2 



The thermodynamic potential ^?qcd(^) is transformed into i^qcol^ + ^nk/S) by the first 
transformation of Eq. (|2]) and it is transformed back into ^2qcd(^) by the Z3 transformation, 
that is, by the second and third equations of Eq. All quantities invariant under the 
extended Z3 transformation, such as the thermodynamic potential and the chiral condensate, 
keep the RW periodicity. Meanwhile, the Polyakov loop <P is transformed as ^ — > ^g-i27rfe/3 
under the transformation (121) and then does not have the RW periodicity. However, this 
problem can be solved by introducing the modified Polyakov loop ^ = <Pe [36] that is 
invariant under the extended Z3 transformation. As an essential property, thus, QCD has 
the extended Z3 symmetry and it is realized as the RW periodicity in the /xj region. 
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Fig. 1: Phase diagram on the 0-T plane predicted by the PNJL model. This diagram has a peri- 
odicity of 9 = 2tt/3. The dashed curves represent crossover deconfinement phase transitions, while 
the solid lines represent the RW phase transition and its Z3 images. Points E are endpoints of the 
RW phase transitions. Point C is a pseudocritical transition point of the crossover deconfinement 
phase transition at 9 = 0. 
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Among many effective models proposed so far the PNJL model is only a realistic effective 



model with both the extended Z3 symmetry and chiral symmetry 36|, |40|, |4l| . As a resu 



of this property, the PNJL model succeeds in reproducing the RW periodicity |36|, 141 1. 
Figured] shows the two-flavor phase diagram in the 6-T plane predicted by the PNJL model; 
the details of the PNJL calculation will be described in Sec. HIl The dashed curves represent 
crossover deconfinement phase transitions, where the pseudocritical temperature at each 9 
is determined by the peak position of the Polyakov-loop susceptibility. 

Roberge and Weiss also showed by using the perturbative and the strong coupling QCD 
that i^Qcol^) is smooth at = 7r/3 (mod 27r/3) when T is low, but not when T is high |42| . 
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This indicates that there exists a phase transition above a temperature Te. This disconti- 
nuity of df2QCD{(^)/d0 is called the RW phase transition. The RW phase transition is known 
to be first order at T > Tg, but the order at the endpoint T = is not clarified yet. The 
presence of the first-order RW phase transition is confirmed by LQCD 4j4l0l|. In Fig. [U 
the solid line represents the RW phase transition predicted by the PNJL calculation. Point 
E is the endpoint of the RW phase transition. Obviously, this phase transition also has 
the RW periodicity. This success of the PNJL model suggests a possibility that the phase 
diagram in the /ir, region is determined from the thermodynamics in the /ii region by using 
the PNJL model, the parameters of which are fitted to reproduce LQCD data in the /ii 
region. Actually an analysis along this line has been made very recently in Ref. |4l|. Thus, 
deep understanding of the phase structure in the /ij region is important to determine the 
phase diagram in the /xr region. 

The study of the /zj region has another important aspect. Roberge and Weiss found that 
at T > Te three Z3 vacua come out alternatively as 6 varies from to 2it and showed that 
at 9 = 71/ 3 (mod 2tt/3) a transition from one of Z3 vacua to another happens 42|]. This 
mechanism is an origin of the RW phase transition appearing at = 7r/3 (mod 27r/3). We 
call this mechanism the RW mechanism. At = 0, one of three vacua is selected. At T < Te, 
there is no Z3 vacua, and hence the RW mechanism does not take place. 

In Fig. [H the dashed curve between points C and E represents a line of crossover decon- 
finement transition. It is natural to think that the behavior of the crossover deconfinement 
transition on the line and particularly at point C is influenced by the critical behavior of 
the RW phase transition at point E. The RW phase transition is studied in Ref. by 
LQCD with four-flavor staggered fermions. The work suggests that the critical behavior 
of the RW phase transition at point E is essential to nonperturbative features of strongly 
coupled quark-gluon plasma (sQGP) appearing at = and Tq < T < 3Tc, where Tq is 
a temperature of point C in Fig. [H Thus, the study on the RW phase transition is also 
important to understand properties of the crossover deconfinement transition and sQGP at 
/i = 0. 

In this paper, using the two-fiavor PNJL model, we investigate properties of the RW 
mechanism and the RW phase transition as a consequence of the mechanism. Concretely, 
the following three points are argued. First, we show that in the RW phase transition the 
charge conjugation (C) symmetry is spontaneously broken, while the extended Z3 symmetry 
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(the RW periodicity) is preserved. Second, we show that the RW phase transition is of second 
order at the endpoint {9,T) = (tt/SjTe). In the /xr region, it is we^ 



1 known that the chiral 

phase transition is of second order at the critical endpoint [13|, |17|, |18|. Singular behaviors 
near the critical endpoints are different between the two second-order phase transitions. 
Third, we argue that the crossover deconfinement transition at /i = is a remnant of the 
RW phase transition at the endpoint. 

In Sec. [ni the PNJL model and the extended Z3 symmetry are explained briefly. In 
Sec. Illlt the RW phase transition and the RW mechanism are analyzed both analytically 
and numerically. Section IIVI gives a summary. 



II. PNJL MODEL 



A. Model setting 



We consider the two-flavor PNJL Lagrangian with fi = iT6, 

C =q{i'j^D'' - mo)q 

+ G,[{qq)' + {qz^,Tqf]-U{<P[A],cP[Ar,T), 



(3) 



where ttiq is the current quark mass, D'^ = d'^ + iA'^ — ijjLdQ and A'^ = SqqA^^ with the gauge 
field y4|^, the Gell-Mann matrix Aa and the gauge coupling g. In the NJL sector, r stands 
for the isospin matrix, and Gs denotes the coupling constant of the scalar-type four-quark 
interaction. The Polyakov potential U, defined in Eq. ([8]), is a function of the Polyakov loop 
^ and its Hermitian conjugate 



1 



ticL, 



-tr,L^ 



(4) 



with 



LU) 



I I (iry44(x, r) 




(5) 



where V is the path ordering, A4 = zAq and = 3. In the PNJL model, <P and <P* are 
treated as classical variables. We also denote <P by Re'^^, where R and (p are the absolute 
value and the phase of ^, respectively. In the chiral limit (mg = 0), the Lagrangian density 
has the exact SU{2)j^ x SU{2)^ x f/(l)^ x SU{3)c symmetry. 
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Using the mean field approximation (MFA), one can obtain the thermodynamic potential 



per unit volume [23|, Q , 



n = -2Nf 



d^p 
(2^ 



+ -ln + 



3i?(p) + -ln + 



(6) 



with 



where Nf = 2, a = {qq) is the chiral condensate, /3 = 1/T and E{p) = a/p^ + with the 
effective quark mass M = mo — 2GsCr. We use W of Refs. |23|, l37[ that has a strong coupling 
inspired form, 

546-"^'^^^* + ln(l + G) 



U = -bT 



with 



G = -6<P<P* + 4(^^ + <P*^) - 3 



*\2 



(8) 



(9) 



The constant parameter a is taken to be 664MeV so as to reproduce the LQCD result on 
pure gauge that the first-order deconfinement phase transition takes place at T = Tq = 



270MeV 



23 



3- 



The vacuum term (the first term of the right-hand side of Eq. ([6])) diverges. It is 
then regularized by the three-dimensional momentum cutoff A: 



P 



1 



(27r) 



— > 



27r2 



dpp'^ 



(10) 



Following Ref. 23|, we regularize the vacuum part, but not the thermal part. Even if the 
thermal part is regularized, it does not change the present result much unless T is much 
larger than T^. 

The parameter b in U means a mixing strength between the chiral and deconfinement 
phase transitions, and is chosen to be 0.015/1^ to reproduce the two-flavor LQCD data in 
which a crossover deconfinement transition occurs around Tq ~ ISOMeV; in ref. 37|], b is 
taken to be 0.03/1^ to reproduce Tq ^ 200MeV in the three-flavor case. 

Hence, the present model has three parameters mo. A, Gg in the NJL sector. Following 
Ref. (igj, we use mo = 5.5 MeV, A = 0.6315 GeV, and G, = 5.498 GeV'^ that reproduce 



the pion decay constant = 93.3MeV and the pion mass = 138MeV. 



In the /ii case, i? and a are real, while is the complex conjugate to <P [36]. Variables, 
X = <P* and 0", satisfy the stationary conditions, 

dQ/dX = 0. (11) 

The thermodynamic potential f2{6) at each 6 is then obtained by inserting the solutions 
X{6) in Eq. ([6]). The thermodynamic potential f2{6) thus obtained does not give the global 
minimum of i7 necessarily. Actually in the case of high T, there exist local minima (unstable 
solutions) in addition to the global minimum (the stable ground-state solution), as shown 
later with numerical calculations. 



B. Extended Z3 symmetry 

The thermodynamic potential f2 of Eq. ([6]) is invariant under the extended Z3 transfor- 
mation, 



* i27rfc/3 



(12) 



This is easily understood by introducing the modified Polyakov loop \P = e*^^ and \P* = 
e~*^<p* invariant under the transformation f|T2|) . The extended Z3 transformation is then 
rewritten as 



(13) 



and i? as 



f] = 12™^ + father ^ 



(14) 



with 



ther 



-2Nf 



2Nf 



(2^3- 



(3 J (27r)^ 



3E(p) + Gsa^, 
ln(l + F) + ln(l + F*) 



U = -bT 



^Ae-^'^m* + ln(l + G) 



(15) 
(16) 
(17) 
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where 

F = SiZ^e-^-^^P) + 3i^*e-2/5'^(p)e3*^ + e-3/5'^(p)e^*^ (18) 

G = -em* + 4(lp'3g-3ie ^ ,^*3g3i0) _ 3(5^5P'*)2_ (19) 

Obviously, i? is extended Z3 invariant, since it is a function of only extended Z3 invariant 
quantities, e^*^ and X = cr, and The explicit 6* dependence appears only through the 
factor e^*^ in Eq. ( IT^ . Hence, if the solutions {X} are uniquely given, they have X = X(e^*^). 
Inserting the solutions back into Eq. f lT4|) . one can see that O = J7(e^*^). Thus, f2{6) and 
X(^) have the RW periodicity, f]{9) = Q{e + 27rA;/3) and X{e) = X{9 + 27rk/3), because 
they are extended Z3 invariant. However, the situation is more complicated at high T, since 
three sets of solutions, {k = 0, ±1), are given; this will be discussed in Sec. IIH B[ 



III. ANALYTIC AND NUMERICAL RESULTS 



A. RW phase transition 

Under the charge conjugation (C), the Polyakov loop and the chemical potential n are 
transformed as — i- <P* and ^ — )■ respectively; for example, see Refs. 2J, |43|]. This 
indicates that the modified Polyakov loop ^ = ^e*^ is also transformed as ip' — )■ ip'*. As 
shown in Eqs. (ITB]) and f[T71) . i7 is invariant under the C transformation. In other words, i? 
is invariant under the transformation 6 — t- —6, if ^ is replaced by This means that the 
solutions X{9) of the stationary conditions ffTTl) satisfy 

ip{e) = ^*{-e), a{e) = a{-d). (20) 

Equation fl20l) indicates that the chiral condensate a, the absolute value \^\ and the real 
part Re[!P'] are ^-even, while the phase ip = arg(!P') = (f) + 6 and the imaginary part Im[!P'] are 
^-odd. Inserting the solutions X{6) into Eq. flT^ . one can see that f2{6) is ^-even. Hence, the 
derivative df2{6)/d6 and the quark number density n = —dO/dfi = —df2/d{iT6) are 6'-odd 
quantities with the RW periodicity. Such 6'-odd quantities 0{9) with the RW periodicity 
satisfy 

lim 0(6 - e) = - lim 0(6 + e) (21) 
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at 9 = 7r/3 (mod 2n/3), because 0(7r/3 - e) = -0(-7r/3 + e) = -0(7r/3 + e). Thus, 0{9) 
is discontinuous when hm^^+o 0(6* + e) is finite. This discontinuity is called the RW phase 
transition, and realized in the high T region, T > = 205 MeV, as shown below. 



-0.02 - 




2 4 6 



e/(7i/3) 

Fig. 2: Thermodynamic potential 17 as a function of 6. The solid curve represents a result of 
the case of T = 400 MeV, and the dashed one corresponds to that of T = 150 MeV. The inset 
represents the ^-dependence of J? at T = 150 MeV in a small scale. 

Figure |2] shows 6'-dependence of il{6) calculated with the PNJL model. At T = 400 MeV 
belonging to the the high-T region T > Tg, the solid curve has a cusp at 6 = n /3 (mod 
27r/3). Thus, df2{6)/d6 is discontinuous there. At T = 150 MeV belonging to the low- 
T region T < Te, meanwhile, the dashed curve is smooth everywhere and has the RW 
periodicity, as shown by the inset. In this low-T case, 6'-dependence of i7 is very weak. In 
the zero-T limit, furthermore, i7 has no ^-dependence, since 

f2\T=o = ^2 = -QNf I ^,E{p) + G.a\ (22) 

where a is obtained by the stationary condition d{Q\T=Q) / da = 0. 

The solution f2{6) is transformed by the charge conjugation C as f2{6) — )■ f2{—6). When 
the solution 0{6) with 6 fixed is considered, C is a symmetry of it only at = and tt; 
note that = tt is identical with 6 = — vr. The ^^-odd quantity 0{6) such as ip and n is 
transformed by C as 0{7i) — )■ — 0(— vr) = — 0(7r) and hence not C-invariant at 9 = n. When 
T < Te, nevertheless, 0{6) is a smooth function of 6, so that it is zero at 6 = it because of 
Eq. (^T^. Thus, we can regard 0{6) as an order parameter of the C symmetry and then use 
ip for this purpose. 

Figure [3] shows the T dependence of ^ at = vr. ip is zero at T < = 205 MeV, while 
it is finite at T > Te. Thus, the C symmetry is spontaneously broken above Te, although 

9 
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Fig. 3: Phase of the modified Polyakov loop as a function of T in the case of ^ = vr. 



it is preserved below Te. Same phase transition takes place at = 7r/3 and bn/S as a 
consequence of the RW periodicity. At T = Te, dip/dT is discontinuous, implying that the 
RW phase transition is of second order there. This result is consistent with the speculation 



of de Forcrand and Philipsen 
made in subsection IIII C[ 



based on LQCD. Further discussion on this point will be 
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Fig. 4: The absohite value of the modified Polyakov loop as a function of T in the case of = vr. 



Figure m shows T-dependence of the modulus \^\ at 9 = n. The derivative d\\l/\/dT is 
also discontinuous at T = Te. This also implies that the RW phase transition is of second 
order at the endpoint {9,T) = (tt/S mod 2tt/3,Te), although is not zero at T < Te and 
then not an exact order parameter of this phase transition. 

Figure [5] presents the modulus as a function of T and 6. The second-order phase 
transition appearing at ^ = 7r/3 becomes crossover as 6 decreases from vr/S to 0. Thus, the 
crossover deconfinement transition at = 0, shown by a rapid change of \^\ with an increase 
of T, is a remnant of the second-order RW phase transition at 6 = n /3. 

The PNJL analyses mentioned above are summarized as the phase diagram in Fig. [H 
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Fig. 5: The absolute value of the modified Polyakov loop as a function of T and 9. 

The temperature of point C (the pseudocritical temperature of the deconfinement phase 
transition at ^ = 0) is Tc = 179 MeV, and temperature of point E (the endpoint temperature 
of the RW phase transition) is Te = 1.1 5Tc = 205 MeV. Thus, Te is higher than Tc. 

B. RW mechanism 

We start with the SU(3) pure gauge system. For ^ = Re^'^, the Polyakov Potential U is 
obtained as 



The potential U has a dependence only through the cubic term cos 30, that is, the 
<p3 _j_ (p*3 ^gj,j^ fjj^g fg^g^ _5 < Q means that U has a minimum at = mod 27r/3, if R is 
not zero. Now we consider the case of = 0. Figure [6] show the R dependence ofU/A'^ at 
= 0. Three cases of T/Tq = 0.99, 1, 1.01 are shown by solid curves from the top to the 
bottom, respectively, where Tq = 270 MeV. For each T, U has two local minima. When 
T < To a global minimum is always located at i? = 0, but at T = Tq the location jumps from 
i? = to 0.45. Thus, a first-order deconfinement transition takes place at T = Tq. Due to the 
i?^ cos 30 term, at T > Tq, there are three Z3 global minima at (-R, 0) = (-Ro^O mod 27r/3), 
where Rq is a value between 0.45 and 1. Hence, the ground state of the pure gauge system 
has a 3- fold degeneracy above Tq, but no degeneracy below Tq. In the pure gauge system. 



U = -hT 54e-"/^i?2 



+ log (1 - + cos 30 - 3i?^) . 



(23) 
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thus, the high-T phase is distinguishable from the low-T one by the number of the vacuum 
degeneracy. 



U/A^ 




Fig. 6: R dependence oiU / A'^ at (j) = ^ mod 27r/3. Three solid curves correspond to the cases of 
T/Tq = 0.99, 1, 1.01, respectively, from the top to the bottom. 

Next, we consider the system with dynamical quarks. As a feature of the system, Q{9) is 
a smooth function of 9 below Te, as shown in Fig. O In order to understand this property, 
we consider the case of small T. First, we assume that ^ tends to zero and a does to a finite 
value (To in the limit of small T. This is a natural assumption and justified below. In the 
case of small but nonzero T, the stationary conditions (fTT]) for <P and a have T dependence 
through factors e"'^^''^^ and e~^". The factor e'^^'-^^ has a maximum e~^^^ at p = 0, and 
near T = the maximum is very close to e = e~^^^° with Mq = itlq — 2G'sO"o. Hence, 
the p integration of 6"^^^^-' is of order e. Further, the factor e~^" is of order because of 
a ~ 2Mq. It is then natural to expand <P and a by powers of e: namely, 'P = (l>ie+<p2^^ ■ ■ ■ and 
a = (Jo + (Tie ■ ■ ■ . One can then derive equations for coefficients <Pj and (Tj of the e expansion, 
order by order, from the stationary conditions ffTTl) . The equations for the leading-order 
solutions, (Tq and ^i, are 

r d^p Mo , , 



= <{ ILL y _r_Z_e-/^^'olP)+/^Mo (25) 



where -E'o(p) = a/p^ + Mq . Equation ( 12^ does not include T and 6, so that the solution 
(Tq is a constant, as expected. In Eq. ( l25l) . tends to zero as T decreases, as expected. 
The absolute value |^i| does not depend on 9, while the phase of (Pi does. Inserting the 
leading-order solutions into the stationary conditions (fTTI) . we can get higher-order solutions; 
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for example, ai = and 



^2 = |2|<?ip + ^ J _|!P_e-2/3i^o(p)+2/3Afo|g2^e^ (26) 



_ 12Nj\<P,\ r d'p Mo pEoipHfSMo (27) 

12iV^ /■ d'p Mo 
Co y (27r)3Eo(p) ^ ^ 



with 



C^.-l-l^'V.G./li^j^-OA (29) 

where ^3 has not been exhibited, because it does not contribute to i7 up to e^, as shown 
below. The second-order solutions, \<p2\ and a2, do not depend on 6, but the phase of 
(p2 and the third-order solution do. Thus, unique set of solutions, a = ^„=o 
<P = Yln=i ^n^"' is obtained at small T. The thermodynamic potential i7 is given by inserting 
the set of solutions into Eq. (j6]). Hence, the potential, i7 = J2n=o ^^"■'e"; thus obtained is 
unique and then smooth in 6; note that \F\ ^ 1 and \G\ ^ 1 and then log(l + F) ^ F 
and log(l + G) ^ G have no singularity. The factors F and G can also be expanded into 
P = E„=2 and G = Y.n=2 Gut"" with 

F2 = 3|<?i|e-^^»(p)+^^^°, G2 = -Q\<Pi\\ (30) 

G3 = -4(3|^i||^>2| -2|<?iHcos(3^),--- . (32) 
Eventually, the f?*^"-* are given as 

12(0) = G^al - 6Nf I ^^o(p), O^'^ = 0, (33) 

12(2) = 2G,aoa2 + 12G3iV^a2 / t^^t^ - ^NfT [ ^^^2 - 6TG2, (34) 

y (27r)3Eo(p) J (27r)3 

^3 ° -4iV^Tcos(3g) / TTTTll^sl 
(27r)'^Eo(p) J {27ry 

-bTGs, (35) 

up to order e^. Thus, i? has no 9 dependence up to order e^, and the ^-dependence appears 
first at order through the factor cos(36'). Therefore, i7 is a smooth function of cos(36') at 
small T, but the 9 dependence is very weak, as explicitly shown in Fig. [2j 
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Above Te, i7 has a cusp at Q = n j?) (mod 27r/3), as shown in Fig. |2l To understand the 
nature of the discontinuity, we consider the case of high T, extend A^/ to be a continuous 
variable and take the small N ^ limit. Since and i?*'^'^'' are proportional to A^j, it can be 
treated perturbatively. To the zeroth order of A^j, the thermodynamic potential is reduced 
to that for the pure gauge theory. In this gauge dominant case, the thermodynamic potential 
has a global minimum at R = Rq satisfying 0.45 < Rq < 1, a^s shown in Fig. [61 Therefore, 
there exit three unperturbed solutions <Pk = |^fc|e*'''''= = R^e^'^'^'^^^ {k = 0, ±1), since U is 
invariant under the Z3 transformation. 

Inserting the zeroth solutions = R^e^'^'^'^^^ {k = 0, ±1) into H of Eq. f[T^ . one can get 
three kinds of thermodynamic potentials f^k{d) {k = 0, ±1). For simplicity, we consider the 
limit of i?o — ^ 1- The thermal fermionic parts f2^i}^'^^{9) are 



other fn ^ ^^f f 



j ^[in [1 + e-'5^^(p)e^^+^2.fe/3] 



where E{p) = a/p^ + {mQ — 2GsO")^- The chiral condensate a in Eq. ( 136|) is determined by 
the stationary condition dfl/da = 0: 

with 

Equations (E?]) and f l55]) show that a is small at high T, independently of Nj, because the 
factor D{p,T,6) tends to zero in the limit of high T. In the numerical calculations for 
Nf = 2, <Pk tends to e*^'^'^/^ and cr becomes small as T increases. In this sense the present 
analysis is consistent with the numerical results. 

The solutions Ok{6) are smooth, but each solution is periodic only with period 27r. Thus, 
each of the solutions does not have the RW periodicity. This does not mean that a set of 
the three solutions does not keep the RW periodicity. In order to understand this clearly, 
we consider the vanishing quark mass limit of M = where the analytic forms of 
are available: 



father ^ [Li4(e^(''-'"'/=^)) + Li4(e-^(''-2-^/3))l , (39) 



14 



where Li^lz) is the polynomial logarithm defined by Li^lz) = Yl'^=i /n^- As shown in Fig. 
[7] (a), the three solutions Q)^'^^{9) are smooth but periodic with period 27r, and fl^^i^{9) are 
obviously Z3 images of fl^'^^{9). As shown in Fig. [7] (b), if the lowest solution is taken at 
each 9, the ground state (GS) thus connected, i7j^g^(6'), is periodic with period 27r/3 but not 
smooth at = 7r/3 (mod 27r/3). As a result of this mechanism, ^2gs(^) = + ^2gs'(0) +W 
is not smooth at ^ = vr/S (mod 27r/3), that is, the RW phase transition appears there. We call 
this mechanism the RW mechanism. In the high T case, thus, the ground state preserves the 
RW periodicity, i.e. the extended Z3 symmetry, although each Q^j^'^'^{9) does not. Eventually, 
the extended Z3 symmetry is held at any temperature. 




0/(71/3) 9/(71/3) 



Fig. 7: Perturbative solutions to the thermodynamic potential in its thermal part: (a) the smooth 
solutions {k = 0,±1) and (b) the ground-state solution f^Q^"^- In (a), the solid, the dotted 

and the dashed curves show the cases of /c = 0, k = +1 and k = —1, respectively. 



The analyses above based on perturbation are essentially the same as Roberge and Weiss 



did with perturbative QCD within its one- loop approximation |42|, |4J]. Actually, the one- 
loop solutions agree with ilf"^^{9) {k = 0, ±1) of Eq. fl39l) . Thus, at high temperature, the 
system has three Z3 vacua, and they appears alternatively as 9 varies from to vr. This is 
the RW mechanism, and when ^ = the system belongs to one of the Z3 vacua. 

The RW mechanism appears also in the full PNJL calculation free from perturbation. 
However, the situation is more complicated as shown below. The stationary conditions f llip 
give three sets of solutions, cr^ and $k = l^fel^*'^'' {k = 0,±1), as expected. Inserting the 
solutions into Eq. (fl^ . we have three solutions i?^. Figure [8] shows ^-dependence of i?^ 
and (pk in the case of T = 400 MeV. Surely, there exist three kinds of solutions. At each 



15 



9, however, at least one of the three disappears; the solution i7o vanishes in the region of 
O.SStt < 9 < 1.467r, and the region is shifted by either 2tt/3 or — 27r/3 for other solutions. 
Whenever exists, 0^ is about 2iTk/3 almost independently of 9, while depends on 9. 
The ground state ^^gs is composed of in region (I) — 7r/3 < 9 < 7r/3 (mod 27r), f2_i in 
region (II) vr/S < 9 < tt (mod 27r), and I2i in region (III) — vr < 9 < — 7r/3 (mod 27r). 




2 4 

e/(7i/3) 




2 4 
e/(7l/3) 



Fig. 8: 9 dependence of solutions to the stationary conditions in the case of T = 400 MeV: (a) 
represents 17^ and (b) does 0^- The solid, the dotted and the dashed curves show the cases of 
k = 0, k = +1 and k = —1, respectively. 

The RW mechanism are analogous to the Dashen phenomenon 4^ in the so-called 0- 
vacuum. Following Witten's analysis [46] on the Dashen phenomenon, we can discuss the 
spontaneous breaking of the C symmetry in the RW mechanism. The C transformation 
changes the sign of 9. Hence, for the thermodynamic potential f2{9) with 9 fixed, C is a 
symmetry of f2{9) only at ^ = or = tt; note that 9 = 7t/3 mod 27r/3 has the same 
property as = tt because of the RW periodicity. If two of the solutions ^2^(9) cross each 
other at 6' = and 9 = it, each solution is C-violating and C interchanges the two solutions 
there. This Witten's argument on C violation can be explicitly confirmed in this case, as 
mentioned below. The charge conjugation C transforms the sign of ip, and hence ip = is 
invariant under C. In addition, ip = n is also invariant under C because ip = n is identical 
with ip = —IT. The solution J7o(^) is C-conserving, since it has tp = at 9 = 0. Meanwhile, 
f2±i{9) are C-violating solutions, because in these solutions ip is neither nor vr at ^ = vr. 
As shown in Fig. [8](a), the C-violating solutions f2±i{9) cross each other at ^ = tt, while the 
C conserving solution fIo{9) has no crossing at = 0. Thus, the C symmetry is conserved 
at ^ = 0, but spontaneously broken aX 9 = -k. This C symmetry breaking appears also at 
9 = ±71 / 3 as a consequence of the RW periodicity. 
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When T < Te, the ^-odd quantity ijj is zero at 6' = and vr, because ijj is a smooth 
function of 6 satisfying Eq. fl2T]) . Hence, the C-symmetry is preserved there. Furthermore, 
ip = at 6 = kn/S with integ resuh of the RW periodicity. As seen in Sec. IIII Al 

6'-dependence of f2{6) is weak, and then ~ 0. Therefore, (p = — 9 —9 at low T. On the 
contrary, when T > Tg, is almost constant in each of regions (I), (II), (HI). The T > 
regime is thus distinguishable from the T < Te one by the ^-dependence of (or tp). 



C. Order of RW phase transition at endpoint 



Susceptibilities Xij of a, R= and ip = arg(ip') can be written as |23l. |32|. |35| 



Xij = {K = cr,i?, V"), 



(40) 



where 



K 



( 'dHi d^n d^n \ 

da^ dcrdR dcrd^) 

d'^Q d'^n d'^Q 

dRda (9_R2 dRd^ 

d^n d'^a 8iH2 I 

\ dipda di>dR ai/>2 / 



(41) 



is a symmetric matrix of curvatures of i7 and Kij is an (i, j) element of curvature matrix K. 
Matrix elements and K^^ are 6'-odd and then zero at T < Te and 9 = ti /?> (mod 27r/3), 
as shown in Eq. fl2Tl) . Hence, the Xij is obtained there as 

1 



where 



Xij = {K^^)ij {i,j = a,R), 



(42) 



(9(t2 dadR 

aRaa 'aW 



(43) 



At the critical point {9,T) = {tt/3 mod 27r/3,TE), thus, the susceptibilities of a and R 
decouple from that of tp, the order parameter of the phase transition. 

Figures ini^a) and (b) present T dependence of susceptibilities, Xa-a, Xrr and Xtpip: Sut 9 = n. 
The susceptibility xv-v ^ divergent peak at T = Te = 0.2048 GeV. This indicates that 
the RW phase transition is of second order at the endpoint {9, T) = {n/S mod 27r/3, Te) and 
also that a 6'-odd quantity such as ip is an order parameter of the phase transition. This 
result is consistent with the LQCD result j4|. 
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0.205 

T[GeV] 



0.21 



Fig. 9: (a) Susceptibilities of the chiral condensate a (dashed curve), the absolute value R of the 
modified Polyakov loop (dotted curve) and the phase of the modified Polyakov loop (solid curve) 
as a function of T at = vr. Xaa is divided by T^, while XRR aiid Xipi^ are multiplied by lOT^ and 
0.2T^, respectively, (b) Susceptibilities of the chiral condensate a (dashed curve) and the absolute 
value R of the modified Polyakov loop (dotted curve) near the critical temperature Te. Xo-o- is 
divided by T^, while XRR a-nd Xipip are multiplied by T^. (c) Off-diagonal elements K^^ (solid 
curve) and Kfi^ (dashed curve) are shown as functions of T. K^j^ and K^^ are divided by T and 
T^, respectively. 

As shown in Eq. (H2|) . the divergence comes from the fact that K^^ = at the endpoint, 
and the susceptibihties of the 6'-even quantities such as Xo-o- and xrr are irrelevant to K^^. 
Therefore, the divergent behavior does not affect Xuu and Xrr- Actually, Xuu and xrr 
have no divergent peak there, as shown in Figs. [9](a) and (b). There is no a priori reason 
that the transition temperature T^- defined at the peak position of Xo-o- coincides with Te at 
9 = tt/2i (mod 27r/3). However, it is shown that the vector-type four-quark interaction and 
the eight-quark interaction make closer to Te [41 |. 

This second-order phase transition is accompanied by the spontaneous breaking of charge 
conjugation (C) symmetry. The C symmetry is a Z2 symmetry, since fl is C-even (6'-even) 
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and then an even function of C-odd (^-odd) quantities such as Im[!P'] and ip. Thus, the Z2 
symmetry is spontaneously broken at T > Te- 

In the /iR region, the second-order chiral phase transition at the critical endpoint is 
known to be accompanied by a spontaneous Z2 symmetry breaking 13|. However, as shown 



in Refs. 
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18| . the Z2 symmetry is not exact such as the C symmetry at the RW endpoint. 



and the fiat direction of the effective potential at the critical endpoint is not the a direction 
but a linear combination of the chiral condensate a, the entropy density s and the quark 
number density n. Therefore, off-diagonal elements of the curvature matrix K do not vanish, 
and therefore susceptibilities of a, s and n diverge simultaneously at the critical endpoint 
where the determinant of K vanishes. 

As seen in Fig. [H^b), Xo-o- and xrr have cusps at Te. Since K^-^ and Kji^ are 6'-odd, as 
mentioned above, they are zero in the region T < Te where the C symmetry is preserved, 
but become finite in the region T > T-^ where the C symmetry is spontaneously broken via 
the second-order phase transition. Eventually, as shown in Fig. [9]^c), K^j^ and Kr.^ have 
cusps at Te. This means that in principle all the susceptibilities have cusps at Te, because 
they are given by the inverse of the curvature matrix K. However, this singular behavior is 
masked by the divergence in x^^. 

The thermodynamic potential fl of Eq. (fT2|) is a function of variables R, ip and a. Taking 
a minimum of i7 in variation of R and a with ip fixed, we can define the potential surface 
as a function of the order parameter ip. Figure [10] shows the potential surface at 
9 = 71. Panels (a)-(c) show the surface in the cases of T/Te = 0.73, T/Te = 1, 1.05 and 
T/Te = 1.95, respectively. For the three cases, surely, is '^2 symmetric under the 

transformation ip — )■ —ip. In the case of T/Te = 0.73, there is a minimum at ip = 0. This 
minimum can be regarded as a ground state i7gr. The C symmetry is not broken in i7gr. In 
the case of T/Te = 1-95, there are two minima aX ip ^ ±7r/3. These correspond to solutions 
f2±i with ip ^ ±27r/3 + 6 ^ =F7r/3 (mod 27r) in Fig. [8l It should be noted that there is no 
minimum at ~ tt. This is consistent with the fact that in Fig. [8] there is no solution with 
~ at = tt; note that (p = ip — 6. At T = Te, the potential surface is fiat around ip = 0, 
indicating that the RW phase transition is of second order at the endpoint {0,T) = (7r,TE). 
The same is true also at ^ = 7r/3 (mod 2tt/3) as a consequence of the RW periodicity. 

Figure [TT] presents the potential surface at yU = 0. Panels (a)-(c) correspond to cases of 
T/Te = 0.73, T/Te = 0.87, 1 and T/Te = 1-95, respectively. Only one solution i7o with 
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Fig. 10: Potential surface as a function of tp at 9 = n: (a) represents the case of T = 150 MeV, 
(b) does two cases of T = 205 and 215 MeV and (c) does the case of T = 400 MeV. In (b), the 
case of T = 205 MeV (215 MeV) is denoted by the sohd (dashed) curve. 



ip = appears above Te, while one solution i7gr with ip = does below T^. This is consistent 
with the fact that in Fig. [8] there is only one solution with = at = 0. At = 0, thus, 
the phase ip is always zero for any T. This property guarantees the C conservation, but 
does not induce any singularity in the T dependence of physical quantities. However, the 
transition of the ground-state structure from the T < Te regime to the T > Te regime makes 
\^\ singular at the endpoint {0,T) = {tt,T-e) of the RW phase transition and then induces a 
rapid change of even at /x = 0. 



IV. SUMMARY 



We have analyzed the Roberge- Weiss (RW) mechanism and the RW phase transition, 
using the PNJL model. Above Te, three Z3 vacua appear alternatively as 6 changes from 
to 27r. As a consequence of this RW mechanism, the extended Z3 symmetry (the RW 
periodicity) is preserved, but C symmetry is broken at 9 = 7r/3 mod 2tt/3. This is the 
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Fig. 11: Potential surface fi{tp) as a function of ip at 9 = 0: (a) represents the case of T = 150 MeV, 
(b) does two cases of T = 179 and 205 MeV and (c) does the case of T = 400 MeV. In (b) the case 
of T = 205 MeV (179 MeV) is denoted by the solid (dashed) curve . 

origin of the RW phase transition. As an order parameter of the phase transition, we can 
select anyone of 6'-odd quantities; a typical one is the phase ip of the modified Polyakov loop. 
The RW phase transition is of second order at the endpoint {0,T) = {tt/3 mod 271 / 3, T^). 

The QCD system has Z3 vacua above Te, but does not below Te- As a consequence 
of the transition, the Polykov-loop has a singular behavior at the endpoint {0,T) = 
(vr/S mod 2tt/3,Te) of the RW phase transition. The singular behavior induces a rapid 
change of at ^ = 0, as presented in Fig. [51 Thus, the crossover deconfinement transition 
at /i = 0, defined by the rapid change of is a remnant of the second-order RW phase 
transition at the endpoint {0,T) = (vr/3 mod 27c/3,T-e). 

Just above Te, one or two of Z3 vacua emerge at each 6 in the PNJL model, while all 
of them appear at each 6 in the RW prediction based on perturbation. Thus, the RW 
mechanism is seen also in the strong-coupling regime which the PNJL model treats, but it is 
somewhat different from that predicted by Roberge and Weiss in the weak-coupling regime. 
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